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4) Procede de visualisation en temps reel d'ecoulements de flu Ides. 

(57) Le procede de P invention consiste a charger dans une 
memoire de travail une texture precalculee, a deplacer sur 
cette texture une fenetre de projection, et a projeter par 
"copy-block" le contenu de la fenetre sur un espace ecran. 
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A 

PROCEDE DE VISUALISATION EM TEMPS REEL 
D'ECOULEMENTS DE FLUIDES 

La presente invention se rapporte a un precede de 
visualisation realiste et anim6 en temps r6el d* 6coulements de 
5 fluides, et & son dispositif de mise en oeuvre. 

Pour visualiser un ecoulement de fluide, en particulier 
un fluide diphasique, dans des canalisations et des reservoirs, 
on pourrait avoir recours au dessin anim& de type cinemato- 
graphique, mais un tel procede ne pourrait etre que preenregistre 
10 et non pas etre la representation en temps r6el d* un phenom&ne 
physique r6el ou mod61is6. 

On pourrait egalement envisager la synthetisation 
informatique des images, mais pour parvenir un bon realisme, en 
temps reel, et pour tenir compte de tous les cas possibles, il 
15 faudrait mettre en oeuvre des moyens informatiques complexes, 
onereux et specif iques. 

La presente invention a pour pour objet un procede de 
visualisation en temps reel d* ecoulements de fluides, en parti- 
culier de fluides diphasiques, ainsi que des phenomenes thermo- 
20 hydrauliques lies & ces fluides, et ce, de fagon la plus realiste 
possible, en faisant appel & des moyens peu on§reux, faciles a 
. mettre en oeuvre (calculateur graphique standaurd) et permettant 
de visualiser distinctement tous les phenomdnes lies & ces 
ecoulements . 

25 Le procede de 1* invention consiste a decouper le 

circuit de fluide en surfaces graphiques elementaires, a pro Jeter 
dans I'espace ecran d'un ecran de visualisation une fenetre 
rectangulaire de la dimension d'une surface graphique §lementaire 
qui se deplace sur line texture preexistante stockee dans une 

30 memoire et correspondant 4 1' aspect sous lequel on desire 
visualiser le fluide. On associe k chacun des pixels d'une 
surface graphique elementaire des "bits de visibilit6", et on 
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projette sur l*6cran les pixels de la texture affect&s de leurs 
bits de visibility. 

Un aspect avantageux de 1' invention consiste & eviter 
les discontinuity dans le mouvement, en utilisant une copie de 
la texture de base dans au moins une direction, la fenetre etant 
deplacee de la texture de base vers sa copie, lorsqu'elle est 
entierement dans la copie 9 et pour suit son mouvement, on la 
ramene dans la texture de base a une position £quivalente. 

L* invention sera mieux comprise a la lecture de la 
description detaillee d*un mode de mise en oeuvre, pris a titre 
d'exemple non limitatif et il lustre par le dessin annex£, sur 
lequel ; 

- la figure 1 est un schema simplifig il lust rant le procede de 
base de 1* invention ; 

- la figure 2 est un schema simpliflS il lust rant 1* utilisation de 
bits de visibility par plan pour determiner un aspect ; 

- la figure 3 est un ensemble de diagrammes illustrant 
1' utilisation de copies de texture pour assurer la continuity du 
mouvement horizontal, selon le proc£d& de 1 * invention ; 

- la figure 4 est un diagramme illustrant la generalisation du 
proced<§ selon la figure 3 aux deux dimensions du plan ; 

-- la figure 5 est un diagramme Illustrant le princlpe d'unicite 
de la projection (effacement et mise 4 jour simultanee) selon le 
procede de 1* invention, et 

- la figure 6 est un diagramme illustrant globalement le 
processus d* animation selon le procedS de 1* invention. 

L* invention est decrite ci-dessous en r6f §rence 4 la 
visualisation sur l*§cran d*un calculateur graphlque standard des 
6coulements fluides diphaslques (melange eau-vapeur) dans un 
circuit de ref roidissement de centrale nuclealre, mais il est 
bien entendu qu'elle n' est pas limitee a une telle application et 
qu'elle peut etre mise en oeuvre pour la visualisation d'etats 
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(ecoulement, convection . ..) de divers fluides ou gaz a une ou 
plusieurs phases et/ou a densites differentes, et ce, dans divers 
autres domaines. 

Les moyens mis en oeuvre peuvent £tre un calculateur 

5 graphique standard k memoire d* ecran, et la m&moire contenant 
les textures peut etre celle du calculateur et/ou une autre 
memoire (par exemple celle d'un autre calculateur). 

Dans le cas de 1* application & un circuit fluide de 
centrale nucleaire, on visualise les phenomenes d*6coulement 

10 eau-vapeur du circuit de ref roidissement de la centrale, les 
phenomenes d* ebullition, d f entrainement de bulles ou de gouttes, 

et en particulier les ecoulements en phase accidentelle. 

L* image est alors pilotee par un logiciel de calcul thermo- 
hydraul ique qui n' entre pas dans le cadre du procede de 

15 1* invention. 

Ainsi que schSmatiquement represents en figure 1, le 
procede de 1* invention part d'une texture 1 precalculee. Cette 
texture peut figurer seulement un liquide, ou un liquide dans 
lequel sont incluses des bulles d* air, ou un gaz, ou un melange 

20 de deux phases (liquide-gaz ...). Dans le cas oH cette texture 
doit representer un liquide seul ou un gaz seul, on y represente 
une alternance de zones claires et de zones sombres af in de 
pouvoir visualiser un ecoulement homogene, ces zones etant des 
bandes sensiblement perpendiculaires & la direction gen§rale 

25 d* ecoulement (c* est-4-dire k la direction de deplacement de la 
"fenetre" decrite ci-dessous). Le contraste entre zones sombres 
et zones claires est choisi au gr6 de l'utilisateur, selon les 
caracteristiques du dispositif de visualisation, et peut meme 
etre un degrade a plusieurs tons. 

30 On a represents en figure 1 l f espace ecran visible 2 

sur lequel est visualise le circuit de ref roidissement d' \m 
rSacteur nucl6aire, en coupe axiale, les parois des conduites, 
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reservoirs et autres circuits pouvant etre representes en gris. 
Pour simplifier le dessin, on n* a represents sur l'espace 2 
qu'une portion 3 de conduite. Cette portion 3 est decoupee en 
"noeuds" graphiques element aires 4.1, 4.2, 4.3, . .., qui sont des 
surfaces graphiques element aires & animer. Pour chacun de ces 
noeuds graphiques, l*effet d* animation est obtenu par les 
projections successives dans 1'espace Scran 2 d'une fenetre 
rectanguiaire 5 de la dimension des noeuds, qui se "deplace" sur 
la texture 1. Les deplacements de la fenetre 5 ne sont pas, bien 
entendu, des deplacements physiques, etant donne que la texture 1 
est dSfinie en memoire par les valeurs des pixels la composant et 
que la fenetre 5 est def inie par un ensemble d* adresses de tous 
les pixels compris dans cette fenStre. 

La projection de la f enStre sur les noeuds successif s 
est effectuSe & chaque mise a Jour de 1 • image affichee sur 
l'espace ecran 2 (cycles de calcul). Cette mise & jour peut par 
exemple etre effectuee au rythme de 24 ou 25 images par seconde, 
comme dans le cas des dessins animes. Cette projection est 
avantageusement une operation du type "copy-block", trSs repandue 
en graphisme informatique, dont les parametres d* execution 
principaux sont la designation de l'origine (sur la texture) et 
de la destination (sur 1* ecran), ainsi que le mode de combinalson 
bit a bit entre les pixels d'origine et de destination. 

Le deplacement de la fen§tre sur la texture suit une 
trajectoire qui est fonction des phenomenes a visualiser. Cette 
trajectoire peut etre parametree dans le temps pour obtenir des 
mouvements rectilignes, curvilignes, alternatifs, indeter- 
ministes ... 

La texture precalculee 1 est stockee dans une memoire 
de travail, qui peut etre la memoire centrale, une partie de la 
memoire ecran, une memoire tampon d' ecran (distincte de la 
memoire Scran), ou tout autre dispositif memoire dependant de 
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1* application particuliere du procede de 1* invention (calculateur 
et logiciel associe). 

Application du procede aux noeuds graphiques de forme 
quelconque : 

5 En general, les noeuds 4 visual is er ne sont pas des 

rectangles, et peuvent avoir une forme quelconque, or 1* operation 
precitee de projection ne peut traiter que des rectangles. Selon 
le procede de 1* invention, on definit dans les tables de couleurs 
associees aux pixels de la texture 4 visualiser des "bits de 

10 visibilite" qui determinent ce qui sera vu on non sur I'ecran, la 
projection s* ef f ectuant sur un rectangle enveloppe du noeud 4 
visualiser. 

On a illustre en figure 2 cet aspect du procede de 
1* invention dans le cas d*une image d'ecoulement de fluide dans 

15 une conduite coudee. On a represents sur la figure 2 plusieurs 
plans PI, P2, P3 . ..Pn dont chacun a une fonction de visuali- 
sation specif ique. Ainsi, le plan PI est celui determinant la 
visibilite des informations relatives £l la parol de la conduite 
coudee. Ce plan PI comporte quatre zones, a savoir deux zones de 

20 parol 6, 8, une zone 7 comprise entre les zones 6 et 8 (interieur 
de la conduite), et deux zones 9, 10 a I'exterieur des parois, 
d'un cote et de 1' autre de la conduite, II suffit done dans ce 
cas de deux bits de visibilite pour donner un aspect particulier 
a chacune de ces zones. Ainsi, par exemple pour les plans PI et 

25 P2, les zones 6 et 8 pour etre visibles en couleur uniforme, 
grise par exemple, sont affectes des bits de visibilite 10, la 
zone 7 des bits 01, les zones 9, 10 des bits 00 Cen noir) et la 
zone de fond d* ecran 11 des bits 11 (en noir egalement). 

Le plan PI est le plan de validation de visibilite des 

30 parois 6, 8 (seules les zones 6 et 8 ont le bit de visibilit6 1), 
le plan P2 celui de validation de visibilite de l'ecoulement 
(seule la zone 7 a le bit de visibilitfe 1), les plans P3 & PN 
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sont les plans de definition de 1* 6coulement, c' est-i-dlre de 
codage des pixels de texture, comprenant les motifs (bulles, 
vagues, ondelettes — ) et le fond. D'autres plans supplemen- 
taires peuvent servir, le cas echeant, 4 d6finir d'autres motifs 
ou elements a visual iser. 

A la visualisation, les plans sont superposes et la 
projection de la fenetre respecte la combinaison de leurs bits 
qui se fait de fagon a obtenir 1* aspect voulu pour les diver ses 
zones. 

Pour pouvoir raettre en oeuvre les techniques de 
projection mentionn§es ci-dessus, il est necessaire que le 
calculateur dispose de tables de couleurs programmables af in de 
pouvoir definir 1* aspect des zones d6termin6es par les bits de 
visibilitfe. 

La complex! te geometrique du noeud graphlque ou la 
richesse de la texture est sans Influence sur la vitesse du 
traitement et la capacity memoire necessaire, du fait que la 
texture du fluide est precalculee et que la forme de ce noeud est 
pr6d6termin§e et memorisee une fois pour toutes. 
Traitement de la fenetre en limite de texture : 

Les dimensions relatives texture/fenetre, alnsi que les 
vltesses d* animation sont a priori quelconques. La situation oH 
la fenetre n'est plus contenue totalement dans la texture se 
presente done souvent. Dans cette situation, 11 faut soit 
admettre des discontinuites dans le mouvement (provoquant une 
impression vi sue lie de "sauts"), soit traiter d'une maniere 
distincte la partie de la fenetre se trouvant dans la text\a-e et 
la partie hors de la texture, ce qui est en g£n6ral incompatible 
avec la frequence de mise a jour de V image necessaire pour 
obtenir une animation reallste. 

La solution selon 1* invention conslste 4 utiliser une 
copie de la texture, adjacente dans la direction du d£placement 
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de la fenetre, 4 la texture origine. Comme represents en figure 

3, la fenetre 11 suit par exemple une trajectoire rectiligne 

horizontale. Dans cette figure 3, les diffSrents instants x n-1 » 

x , x „ , correspondent a des instants d* iteration successif s du 
n n+l 

calculateur commandant la visualisation. En haut de la figure, a 

1* instant x « , la fenetre 11 commence 4 sortir par la droite de 
n- 1 

la texture origine 12 et sa partie droite se trouve en fait dans 
la texture 13, copie de la texture 12, qui est disposee 4 la 
droite de la texture 12, adjacente 4 celle-ci. A 1* instant x^, la 
fenetre 11 est presque entierement contenue dans la texture 13 
(milieu de la figure 3). Juste avant 1* instant * n+1 » la fenetre 
11 a fini de quitter la texture origine 12, et 4 1* instant * n+1 
(bas de la figure 3), la fen§tre 11 est ramenSe dans la texture 
origine 12 4 une position equivalente 4 celle qu'elle aurait dans 
la texture copie 13, et ainsi de suite. On obtient ainsi une 
continuity dans 1* eff et visuel, et un traitement du contenu de la 
fenetre tou jours similaire. 

Le calcul, en fonction du temps, de la coordonnSe 
(i=n-l, n, n+l, ...) pointant 1* origine de la projection de la 
fenetre sur la texture est particulidrement simple : c* est le 
reste de la division entiere du chemin total f ictif parcouru par 
la fenStre, par la dimension, dans le sens du mouvement, de la 
texture, solt : 



x^ = modulo 



k * numero de cycle i * vitesse 4 reprSsenter 
dimension de base de la texture 



expression dans laquelle k est un coefficient d* adaptation 
proportionnel 4 1' increment de temps r6el par cycle de calcul, et 
i, le numSro de cycle de calcul, qui est increments 4 chaque 
interatlon d* image . Tous les termes de 1* expression sont entiers, 
et x^ peut §tre directement une coordonn§e de pixel de 1* image 4 
reprSsenter. 
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On a represents en figure 4 la generalisation du 
proceed de 1* invention & une trajectoire quelconque. La texture 
origine 14 est copiee quatre fois : une copie 15 sur sa droite, 
et, en-dessous des textures 14 et 15, les textures 16, 17 sont 
respectivement les copies des textures 14 et 15. Ainsi, on 
obtient la continuity d 9 animation aussi blen pour une trajectoire 
(Tl) sensiblement horizontale, que pour une trajectoire (T2) 
sensiblement verticale ou pour une trajectoire (T3) oblique. Bien 
entendu, pour un mouvement en sens inverse sur une trajectoire 
slmllaire k Tl, T2 ou T3, les roles des textures sont inverses 
par exemple pour Tl parcouru en sens Inverse, la texture origine 
serait 15 et la copie serait 14. 

Unicite de 1* operation de projection dans l*espace ecran : 

Lors d'une iteration pour tin noeud de 1'espace Scran, 
il faut rSaliser les deux fonctions sulvantes : 

1. effacement de la texture projetSe dans le noeud au cycle de 
calcul precedent ; 

2. ecriture dans le meme noeud de la nouvelle texture du cycle 
courant . 

Ces deux fonctions doivent etre real i sees 
simultanement, par une unique operation de projection, pour 
' eviter les effets visuels de cllgnotement des motifs de la 
texture et pour benfef icier d'une frequence d* animation e levee. 

Selon 1' invention, ces conditions sont remplies grace k 
la memorisation, dans la memo ire de travail, de deux textures de 
motifs identiques, mais de codages complement aires, les motifs 
etant codes, dans le cas present, sur deux bits pour chaque 
pixel. De plus, un op6rateur de projection different est associe 
a chacune de ces textures. Par consequent, pour un motif donnS, 
la memo ire de travail contient deux textures, par exemple : 
- texture normale : codage motif : 11, operateur de projection : 
OU 



it 
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— texture complementaire : codage motif 00, operateur de 
projection : ET 

Le processus d* animation utilise alternativement 
chacune de ces deux textures avec 1' operateur associS. Ainsi, 
5 apres initialisation du noeud (projection de la texture de base, 
seule), on obtient la sequence suivante. 

Cycle de calcul Texture utilisee Type de projection 

et resultante 

0( Initialisation) normale neant 



10 n-1 complementaire ET 

n normale OU 

n+1 complementaire ET 

n+2 normale OU 
... ••• 

15 On a represents en figure 5, de fa^on tres simplifiee 



un exemple de mise en oeuvre de cet aspect de 1* invention. Sur 
cette figure, on n'a represents que les deux plans des pixels 
concernant les motifs de la texture* Pour la texture normale, 
representee en haut a gauche, le plan 18 relatlf au motif ne 

20 comporte, pour simplifier les explications, que quelques motifs 
20 dissemines dans ce plan, et le plan 19 est unlformSment a 1. 
Le codage des pixels a visuallser est, 01 pour le fond texture et 
11 pour les motifs. 

Pour la texture complementaire, le plan 22 contient les 

25 roemes motifs 23 que ceux du plan 18 et le plan 21 est 
uniformSment a 0. 

Ainsi pour la texture normale, 1' operation 0U 
appliquSe, sur une texture complementaire sur l'ecran, donnera 
"ll" pour les motifs et 01 pour le fond c*est-a,-dire restltue la 

30 texture normale. 

Inversement, pour la texture complementaire, 
I' operation ET appliquSe sur une texture normale sur 1' Scran 
donne 00 pour les motifs et 01 pour le fond, c'est-a-dire 
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restltue la texture compliment aire. 

AInsi, sur le diagramrae de la figure 5, si, pour une 
initialisation du premier cycle CI, le calculateur commande 
1* operation "ET", la texture complementaire est visible stir 
1* ecran. Au cycle suivant C2, apres deplacement de la fenetre, 
l'operateur 0U est utilise 4 la projection. La texture normale 
vient remplacer sur 1* ecran la texture complementaire dans le 
noeud avec deplacement relatif dans la direction de deplacement 
des motifs de texture par rapport A la position precedente. Pour 
le cycle suivant C3, l'operateur de projection ET fait apparaitre 
sur 1'espace ecran la texture complementaire qui vient remplacer 
la texture normale avec un deplacement elementaire des motifs 
dans la direction de deplacement, par rapport 4 ce qui etait 
visible sur 1'espace ecran au cycle C2, et ainsi de suite. Si 
1'alternance des projections, del imi tees par la fenetre de 
projection, des textures normale et complementaire se fait 5. un 
rythme suffisamment rapide (par exemple toutes les 80 ms) t et si 
le pas des mouvements relatifs de la fen§tre de projection est 
suffisamment fin (de l'ordre de quelques pixels), on obtient un 
effet d* animation realiste et non saccade. L* utilisation des 
textures normale et complementaire avec projection ET/OU alternee 
permet de regrouper en une seule operation de calcul les 
fonctions d'effacement et d'ecriture des motifs du noeud 
correspondant . 

Les bits des plans de textures qui ne concernent pas la 
definition des motifs k deplacer ont une valeur telle que les 
plans correspondants du noeud dans 1'espace ecran ne sont pas 
affect§s par les projections d' animation, 4 savoir "0" pour la 
texture normale (operation 0U) et tt l" pour la texture complemen- 
taire ( opera teur ET). Les bits concemes dans ces plans de 
textures sont les bits de visibility et les bits de codage de la 
couleur de fond de texture. 
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Sur le diagramme de la figure 6, on a represents de 
maniere tres schematisee le processus d' animation d'un noeud de 
1 9 image k visualiser. 

Sur la partie supSrieure du dessin (se rapportant & une 

5 partie de la memo Ire de travail), on a represents une texture de 
base normale 18 et une texture de base complSmentaire 19 (on a 
represents en traits interrompus les dSbuts des copies de chacune 
de ces textures de base). Sur la partie infer ieure de la figure 
6, on a represente une partie de l'e space Scran cor respondant , et 

10 en particulier un noeud 20 d'une portion de conduite coudSe 21. 

Dans la texture 18, on a schema tiquement represents 
quatre positions successives 22A & 22D d'une fenStre se dSpla$ant 
le long d*une trajectoire 23 sur cette texture. Les positions 22A 
& 22D sont representees SloignSes les unes des autres pour la 

15 clartS du dessin, mais en rSalitS elles sont presque confondues, 
le dSplacement d'une position & la suivante pouvant correspondre 
a un pixel sur l'espace Scran. Les positions correspondantes sur 
la texture 19 sont rSfSrencees 24A & 24D, sur la trajectoire 25 
(identique d. la trajectoire 23). 

20 On a reprSsentS en traits continus les positions dans 

la texture qui correspondent & une opSration de projection de la 
fenStre, a savoir 22A et 22C dans la texture 18, et 24B, 24D dans 
la texture 19. Les autres positions, sont reprSsentSes en traits 
interrompus et correspondent a une projection depuis 1* autre 

25 texture que celle considerSe. 

Les projections visibles sur l'espace Scran sont 
rSfSrencSes VI 4 V4 et sont faites dans l'ordre VI & V4. Pour VI 
(opSrateur OU pour 18), la fenetre 22A est projetSe sur le noeud 
20. Pour V2 (operateur ET pour 19) la fenetre 24B est projetSe 

30 sur le noeud 20, et ainsi de suite pour les positions suivantes 
des fenetres dans 18 et 19. Si les positions successives des 
fenStres sont tres proches les unes des autres et si les 
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operations de projection sont commandees a un rythme suff isamment 
rapide (par exemple au mo ins 15 fois par seconde), on obtient une 
animation realiste et non saccadee. 
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RE VEND I CAT I ONS 

1. Proced§ de visualisation en temps r£el d' ecoulements 
de fluides et de phenomdnes thermohydrauliques 116s k ces 
fluides, caractferise par le fait qu* il consiste k decouper le 
circuit de fluide en surfaces graphiques Element aires, k pro Jeter 
5 dans l'espace ecran d*un Scran de visualisation une fen€tre 
rectangulaire de la dimension d'une surface graphique 616mentaire 
qui se deplace sur une texture preexistante stockee dans une 
memo ire et correspondant k V aspect sous lequel on desire 
visualiser le fluide. 

10 2. ProcedS selon la revendication 1, caract6ris6 par le 

fait que l'on definit dans les tables de couleurs associees aux 
pixels de la texture k visualiser des bits de visibility qui 
determinent ce qui sera vu ou non sur 1* ecran, dans le but de 
conserver 1* image de la surface graphique 61ementaire lors de 

15 1* operation de projection. 

3. Proced§ selon I'une des revendications 1 ou 
2, caracteris6 par le fait que l f on utilise une copie de la 
- texture de base dans au moins une direction, la f enetre £tant 
deplac&e de la texture de base vers sa copie, lorsqu'elle est 
.20 entidrement dans la copie et poursuit son mouvement, on la ramdne 
dans la structure de base a une position equivalente. 

4. Procede selon l'une des revendications precedentes, 
caracteris§ par le fait que l*on memorise dans la m6moire de 
travail deux textures de motifs identiques, de codages 
25 compl&nentaires, chacune des deux textures servant 
alternativement k la projection sur l'espace ecran. 
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5. Proced6 selon la revendication 4, caracteris6 par le 
fait que l*on assocle un operateur de projection different & 
chacune des textures. 

6. Procede selon la revendication 5, caracterise par le 
fait que l'un des operateurs est un "ET" logique, et 1* autre un 
n 0U" logique. 
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